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Abstract 


We revisit the calculation of the three-loop diagrams for the radiative neutrino mass generation 
and consider some relevant constraints on the model recently proposed by L. Jin, et al (2015) |l|. 
We find that the previous approximation is inappropriate due to the neglect of some important 
contributions, and the benchmark point proposed can neither give rise to enough neutrino masses 
nor accommodate these additional constraints, such as the validity of the perturbation theory, the 
electroweak precision measurements, and the neutrinoless double beta decays. 
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The smallness of neutrino masses and the existence of dark matter are the two phenomena 
that require physics beyond the standard model. It is more intriguing that these two physics 


have a common origin. In Ref, 


■ Q. 


Jin, Tang and Zhang proposed an interesting model which 
generates neutrino masses at three-loop level, trying to explicitly provide such a connection. 
However, in the original estimation of neutrino masses, some important contributions, i.e. 

were neglected. Here, we revisit the 


Feynman diagrams (b) and (c) in Fig. 1 of Ref. [l|, 
calculation of neutrino masses numerically by including these additional contributions, and 
hnd that their effect on the hnal results is so large that the approximation used previously 
cannot be justihed. Then, we consider some further constraints on the model, such as 
the validity of the perturbation theory, the electroweak precision tests (EWP'! 
neutrinoless double beta (0z//5/5) decays, which have not been discussed in Ref. 


s) and the 


IK- 

By explicit numerical calculations of the three kinds of Feynman diagrams in Fig. 1 of 
Ref. [l| with the given benchmark parameters, we obtain the following neutrino mass matrix 


rrii, = 


0.97 1.97 0.68 
1.97 6.20 5.43 
0.68 5.43 7.91 




X 10 


-12 


GeV, 


( 1 ) 




which is already excluded by the neutrino oscillation data, e.g., quantitatively {my)rT should 
not be smaller than 10“^^ GeV at 3 ct G.L. Note that our calculations include the contribu¬ 
tions from the diagrams (b) and (c), which were previously neglected in Ref. [l|. The rational 
behind this approximation is that the mass of the neutral scalar A°, = 4331 GeV, for 

the choosing benchmark point is so heavy that it greatly suppresses the diagrams (b) and (c). 
However, our numerical investigation clearly shows that this argument is not appropriate. 
In order to quantify the effects of diagram (b) and (c), we specify the contributions to the ee 
element from three diagrams = (2.25, — 0.35, — 0.93) x 10“^^ GeV, respectively. 

Obviously, the sum of the diagrams (b) and (c) gives the same order of magnitude as the 
diagram (a), but with an opposite sign, which leads to a large cancellation as well as the 
smallness of the neutrino masses for the benchmark point. Thus, it is not valid to only 
consider diagram (a) when estimating the neutrino masses for the present model, as was 
done in the original paper. 

Now we turn to some additional constraints on the present model which were previously 
omitted. Firstly, note that the mass difference between the two charged scalars, H^ 2 i can 
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be determined to be 


= 


\-jV^! S2/3 , 


( 2 ) 


where S 2/3 = sin(2/3). For the benchmark scalar masses and mixing angle, we can obtain the 
absolute value of the dimensionless coupling constant IA 7 I = 357, which is too large to be 
allowed by the perturbation theory [^. 

Moreover, EWPTs have already placed strong constraints to the present model. In par¬ 
ticular, the most stringent one is given by AT G [—0.04, 0.12] at 95% C.L. j^, while the 
extra scalars give new contributions to T as follows. 



( 3 ) 


where 



( 4 ) 


With the above formula, the benchmark point leads to AT = —116, which is clearly outside 


the experimentally allowed range. 

It is also interesting to note that this model can induce a large 0i'f3f3 decay rate via 
the so-called short-distance channel shown in Fig. [H which is expected to be much larger 


than the conventional long-distance one In the present model, such a short-distance 



FIG. 1. On/3/3 from short-distance channel in the present model 


contribution to the half lifetime for Oz//3/3 decay is given by 




( 5 ) 


where rup is the mass of the proton, Gqi the phase space factor. 



( 6 ) 
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with 


n = 


r = - 

a 

2 

IGtt^ 


r(3) f' 
Jo 
r ^ 

/ 

i=i 


'1—n 


dx^ 


dxo 


X1X2 


+ X2m]j^ + (1 - xi - X2)m|rJ^ ’ 

X1X2 

[xim]j^ + X2m]j^ + X2.m%. + {l- Xi- X 2 - X3)m^o]^ ’ 


(7) 


and A ^3 the nuclear matrix element enve 
numerical values of Gqi and in Refs. 


OTing the operator UL'y^diUL'jfidLeRe^j^. Given the 
LZI, we can predict the expected half lifetimes for 
various conventional targets. Unfortunately, with the present benchmark, the predicted half 
lifetimes for all the target nuclei are already excluded by the current experimental limits. 
The largest discrepancy comes from the measurement of ^^®Xe, with the lower bound on the 
half life of 1.9 X 10^® yr [4, LZI], which is compared with the calculated one of 1.1 x 10^® yr. 

Finally, an interesting question is whether we can find some new benchmark point for 
the present model satisfying all the constraints. Note the relation “ 

If we restrict to A 7 < 5 allowed by perturbativity arguments and consider the 
case with only two Z 2 -odd Majorana fermions, we hnd that the {mu)ee obtained is always 
smaller than 10“^^ GeV, which is already excluded by the neutrino oscillation data with the 
vanishing lightest neutrino mass. From this viewpoint, we conclude that it is difficult in 
hnding a viable benchmark point. 
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